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Abstract 

Recently, the visibility graph has been introduced as a novel view for analyzing time series, 
which maps it to a complex network. In this paper, we introduce a new algorithm of visibility, 
" cross- visibility" , which reveals the conjugation of two coupled time series. The correspondence 
between the two time series is mapped to a network, "the cross- visibility graph", to demonstrate 
the correlation between them. We applied the algorithm to several correlated and uncorrelated 
time series, generated by the linear stationary ARFIMA process. The results demonstrate that 
the cross-visibility graph associated with correlated time series with power-law auto-correlation 
is scale-free. If the time series are uncorrelated, the degree distribution of their cross-visibility 
network deviates from power-law. For more clarifying the process, we applied the algorithm to 
real-world data from the financial trades of two companies, and observed significant small-scale 
coupling in their dynamics. 
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I. INTRODUCTION 



Coupled series, and series of long range correlations could not be thoroughly understood, 
unless a global view of them is represented Complex networks {2^4], on the other hand, 
have provided a global understanding of multi-component systems. Recently, the idea of 
complex networks has been implemented for single variable systems , and mapping 
between time series and networks is proposed. For example, nodes in such networks can be 
extracted by binning the time axis or the level axis 

While some features of time series are not simply observed with standard approaches, they 



can be easily studied when converted to a network [8( . For example, financial trades, which 
their dynamics seem to be alike to a white noise, include memory |9l-|ll| and represent high 
clustering when converted to complex networks. Moreover, investigating a map between time 
series and networks is a method of constructing several prototypes of complex networks 12]. 

According to a recent work [3], the visibility algorithm has been introduced as a fast 
computational algorithm of mapping a time series to a network. Several properties of the 
time series are served in the resulting network's structure. For example, a periodic series 
results in a regular lattice, a fractal series results in a scale-free network and a random series 
is in correspondence with a random graph. 

Visibility is the representation of the question that how much local information is attained 
from a global structure. Based on this approach, new global properties are introduced for 
time series. For example, along with scale-free properties, small-world characteristics [y| 



13] . Also, in order to 



15] is proposed. Vis- 



are demonstrated for the fractional Brownian motion time series 
distinguish randomness in time series, a modified visibility algorithm 
ibility algorithm is found to be an appropriate approach in different interdisciplinary areas. 
For example, the method is pertinent to social studies [16!]. In addition, mathematical struc- 
tures such as fractal and complex chaotic time series have been studied using the visibility 
algorithm 

In this work, we introduce a new visibility algorithm, based on the mutual information 
of two time series. We call this type of visibility the cross-visibility. According to this 
method, the visibility criterion is represented for two time series, i.e., a time series looks at 
its components through the obstacles of a second time series. This approach is especially 
appropriate for analyzing real-world data with power-law autocorrelation. The algorithm 
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provides a framework to understand how complex time series lead each other, and deter- 
mines the direction of information flow between them. In addition, coupling in different 
scales is demonstrated. Investigating the cross-visibility between well-known correlated and 
uncorrelated time series could be a criterion for characterizing the cross-visibility. For this 
purpose, we studied the fractional Gaussian noise (fGn) time series. 

Previously, it has been demonstrated that the Hurst exponent for the fractional Brownian 
motion series is distinguished by the structure of the degree distribution in their visibility 
graph We extended this analysis for correlated fGn time series, and characterized the 
power-law exponent in the resulting cross-visibility networks' degree distribution, according 
to their Hurst exponent. Also, in order to depict the implementation of the approach to real- 
world data, we used the data from recent financial trades for two companies. We depicted 
the information flow in several scales between pertinent data. 



II. THE CROSS- VISIBILITY ALGORITHM 



According to the standard visibility algorithm [13j], every component in the time series 
{yi}, with i = 1 ... N, is mapped to a node of a graph. Node i and node j {i < j) are 
connected, if they are visible to each other without any obstacles between them: 

i<Vk<j; y k <y l + ^- Z ^(k- i ). (1) 

J -« 

If node i is visible to node j, node j will also visible to node i. Therefore, the corresponding 
network is non-directional. Elements in the time series with amount high above the others 
have high degrees, and create the singular nodes of the network. 

The cross-visibility algorithm maps two time series, {x{\ and {yi}, into two different 
networks, the cross-visibility networks. The algorithm consists of two major steps: 

a. In order to make the time series comparable, they are normalized to reveal dimen- 
sionless variables. If the sequences are stationary, they are normalized to their mean and 
variances. Therefore, new sequences are generated as {£i} = {^-^} and {yi} = {^-^}. 
x and y are the mean values of the series {xi} and {yi}, respectively; a x and a y are the 
corresponding variances. If the series are non-stationary, they are normalized to their cor- 
responding maximum value. 
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b. Every component in the time series {x{\ is mapped to a node of a graph. Node % is 
connected to node j, if: 

i<Vk<j; yk^Vi + ^^ik-i) (2) 

or 

i<\lk<j- y k >y l + ^2^1(k-l). (3) 

This rule could be interpreted as node i looking at the components of {%} time series, 
through the obstacles of the shifted time series {yj} = {yj — + x^}. The structure of this 
network demonstrates how {xi} leads the time series {?/«}. Equation prefers to the visibility 
from the top view, and equation [3] demonstrates the visibility from the beneath view. We 
depict this graph as G xy . However, another network, G yx , could be generated by changing 
the role of {x^ and {?/,}• 

The cross-visibility can bring a directional measure of mutual information between the 
time series. If the time series {yj} is permuted randomly to construct a time series ( {y^}) 
uncorrelated to {xj}, the structure of the cross- visibility network G x,yP deviates from the 
original one. Each network could be evaluated through its degree distribution. 



A. Characterizing cross-visibility 

In order to interpret the results from the cross- visibility algorithm, it is worthy to analyze 
the time series that are confirmed to be coupled or uncoupled. We investigated the time 



19H21j process. The analysis is a criteria 



series generated by the linear stationary ARFIMA 
for further investigations. 

Through a stochastic recursive rule, the ARFIMA process generates a power law auto- 
correlated series {yu}'- 

oo 

yi = ^2 a j(p)yi-j + (4) 

where, a,j(.) is a weight, defined due to an independent real variable < p < | as a,(p) = 
r(-p)r(i+j) ■ Vii the error term, is a member of an i.i.d sequence, which consists of random 



normal distribution samples. The Hurst exponent, H, is related to p as H = 0.5 + p 20 ]. 



These time series depict power-law autocorrelation 22j . and If exactly same error sequence 
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FIG. 1: (Color on-line) Degree distribution for several cross-visibility networks corresponding to 
different Hurst exponents is depicted. (Blue) squares demonstrate the result for the coupled series, 
(red) stars are for the uncoupled series, and (black) triangles are for the permuted series. Coupling 
between several series with H E {0.6,0.7,0.8,0.9,1} is investigated. 0.6 — 0.7 and 0.7 — 0.6 are 
included in a and b, 0.6 — 0.8 and 0.8 — 0.6 are included in c and d, 0.6 — 0.9 and 0.9 — 0.6 are 
included in e and f, and 0.6 — 1 and 1 — 0.6 are included in g and h, respectively. Both axes are 
sketched in logarithmic scale. 

is used in generating the two time series, they show power-law cross correlation. In contrast, 
independent error terms reveal uncoupled time series. 

We applied the cross-visibility algorithm to coupled time series, with different indepen- 
dent terms p G [0,0. 5[, which are in correspondence with Hurst exponents H G [0.5, 1[. 
FIG. [U depicts the degree distribution associated with each network. A power- law degree 
distribution is demonstrated for the coupled series. The cross-visibility network for the 
uncoupled time series have completely different structure in contrast to the corresponding 
network of coupled series (stars and triangles in FIG. [T]). Moreover, the average degree of 
the distributions is higher for the uncoupled series. 

If the time series are correlated, the corresponding exponent of the power-law degree 
distributions, depicted as j xy for G xy , determine how the series are leading each other. ^ xy 
for several series with Hurst exponents H G [0.5, 1[ is depicted as a color plot in FIG. [2j In 
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FIG. 2: (Color on-line) The color plot demonstrates the exponent of the power-law degree distri- 
bution associated with the cross- visibility networks of coupled series with H £ [0.5, 1[. The result 
is sketched for the cross- visibility networks G Hl ' H2 , where Hi and H2 are the corresponding Hurst 
exponent for the first and the second time series, respectively. 
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FIG. 3: (Color on-line) Power-law exponent for the degree distribution of cross-visibility graphs 
associated with coupled time series is sketched as a function of the difference between the cor- 
responding Hurst exponents. The result is depicted for different Hurst exponents, Hi, for the 
cross- visibility network G Hl,H2 . The power-law exponent decays linearly with H2 — H\. 
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FIG. 4: (Color on-line) The closing price of simultaneously recorded data for Microsoft and IBM 
is analyzed by the cross-visibility algorithm. Degree distribution associated with the original data 
is depicted by (blue) squares. (Red) dots depict the degree distribution for the cross- visibility 
network of permuted time series, G x,yP . The deviation is mostly concerned with the small scales, 
and also the pick of each distribution is shifted for uncoupled time series. Both axes are sketched 
in logarithmic scale. 

order to determine the exponent of the degree distribution we used a maximum likelihood 

n 

estimation (MLE) [23[. In FIG. [3] a linear dependency between jH^fo and H 2 — Hi is 
depicted. This result demonstrates that the Hurst exponents could be distinguished by 
evaluating the power-law exponents. The result is sketched for several Hurst exponents Hi. 



B. Application in empirical data 

For real-world data, it is important to confirm if the time series are coupled and to 
what extent they follow each other. In order to determine the coupling between the two 
time series, {x{\ and {yi}, we put the associated degree distribution of their cross-visibility 
networks into contrast with the corresponding network of permuted time series. 

We analyzed the cross-visibility between the closing price of two companies, Microsoft 
and IBM. The cross-visibility method is applied to the normalized logarithm of the prices. 
For each company, return prices are stationary and are normalized to their average and 
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FIG. 5: (Color on-line) General Hurst exponents IBM, Microsoft and their coupled estimated using 
multi- fractal detrended cross-correlation analysis. 

standard deviations. The degree distribution of the cross-visibility graphs associated with 
original series and permuted time series is depicted in Fig. HJ For the original series, the 
associated degree distribution is a power-law. However, if one of the series (the second series) 
is permuted, the associated degree distribution deviates from the original distribution. The 
most significant deviations are concerned with the small degrees. This however means that 
there is strong correlation in small fluctuations. Deviation in the tail of the distributions 
corresponding to high degrees is negligible. These time series do not follow each other in 
their high amplitude fluctuations. The value of the degree corresponding to the pick of the 
distribution is a characteristic of the lag time response between the two time series. 

It is worthwhile to compare the results with the systematic obtaining of the cross- 



correlation analysis (DXA) 



24j-|26j|. Fig. [5] shows the results for DXA, and indicates how 
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Microsoft leads coupling in higher moments and IBM in lower moments. One important 
difference is that DXA is a symmetric method to analyze the cross-correlation, but the 
cross-visibility introduces two networks of coupling between the two time series. 

III. CONCLUSIONS 

- Understanding the coupling between two time series and how information transfers 
between them is a crucial concept in various areas of complex systems. For this purpose, in 
this paper we depicted this coupling as a network map. 

- In general, as the coupling between two time series is asymmetric, there are two networks 
to present the coupling between them. 

- Our results demonstrate that in average, the number of degrees corresponding to the 
coupled series is less than that of the uncoupled data. There is a shift in P(k) to higher 
degrees in uncoupled time series, which could be a criteria to measure the value of coupling. 

- Our findings show that, increasing the difference between the Hurst exponent of two 
fGn time series, the corresponding degree distribution decays with slower slope in higher 
degrees, which means higher coupling or higher correspondence between the two time series. 
This approach is appropriate for investigating coupling for different empirical data. In 
addition, through this approach the conjugation in different scales between time series could 
be studied. In the end, using this approach, the coupling between Microsoft and IBM has 
been viewed in two coupled networks. 
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